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In this article, we have given a systematic formulation of the new generalized competing mecha- 
nism: the Glauber-type single-spin transition mechanism, with probability p, simulates the contact 
of the system with the heat bath, and the Kawasaki-type spin-pair redistribution mechanism, with 
probability 1 — p, simulates an external energy flux. These two mechanisms are natural general- 
izations of Glauber's single-spin flipping mechanism and Kawasaki's spin-pair exchange mechanism 
respectively. On the one hand, the new mechanism is in principle applicable to arbitrary systems, 
while on the other hand, our formulation is able to contain a mechanism that just directly combines 
single-spin flipping and spin-pair exchange in their original form. Compared with the conventional 
mechanism, the new mechanism does not assume the simplified version and leads to greater influ- 
ence of temperature. The fact, order for lower temperature and disorder for higher temperature, 
will be universally true. In order to exemplify this difference, we applied the mechanism to ID 
Ising model and obtained analytical results. We also applied this mechanism to kinetic Gaussian 
model and found that, above the critical point there will be only paramagnetic phase, while below 
the critical point, the self-organization as a result of the energy flux will lead the system to an 
interesting heterophase, instead of the initially guessed antiferromagnetic phase. We studied this 
process in details. 

PACS number(s): 05.50.+q, 05.70.Ln, 64.60. Cn, 64.60.Ht 



I. INTRODUCTION 



In recent years, there have been continuing efforts towards a clear picture of the self-organization phenomena in 
the phase transitions of magnetic systems. Most of the works Jl]-§| have been concentrated on Ising and Ising-like 
systems, governed by two competing dynamics: Glauber's single-spin flipping mechanism Q and Kawasaki's spin-pair 
exchange mechanism [To| l, both with a probability. The system is coupled to a heat bath at a given temperature 
while being subject to an external energy flux. Glauber's flipping mechanism is to simulate the contact of the system 
with the heat bath. Changing the order parameter, it favors lower system energy. On the other hand, Kawasaki's 
exchange mechanism is to simulate the influence of the energy flux. Keeping the order parameter conserved, it favors 
higher system energy. With these two competing mechanisms and the corresponding master equation, people expect 
to obtain the evolution of the system. As exact treatment is not possible for 2D Ising model, consequently Monte 
Carlo simulation and methods as Dynamic Pair Approximation have been employed. The results obtained helped 
to determine the interesting phase diagrams. (However, people are surprised to find there contradictions in the 
predictions of MC simulations and the dynamical mean-field theory, since they are both proved successful to yield 
good qualitative results in other studies. Though revisions of MC simulations are made and more accurate versions 
of mean-field theory are presented, the puzzle still remains ||.) 

In our earlier studies, we have presented single-spin transition mechanism [ fil]Jl2]| and spin-pair redistribution mecha- 
nism Jl3| . These two dynamics are natural generalizations of Glauber's single-spin flipping mechanism and Kawasaki's 
spin-pair exchange mechanism respectively. They have similar mathematical expressions, and become counterparts of 
each other in the non-conserved and conserved processes respectively. As an example of the applications, we studied 
kinetic Gaussian model with both of them respectively. Our study shows that, in translational-invariant lattices, 
the dynamic critical exponent z = 1/v = 2 is independent of space dimensionality and the governing dynamical 



mechanism. Its dynamic properties are summarized in Sec. [V 



In this article, we formulate the competing dynamics combining single-spin transition and spin-pair redistribution. 
As these two mechanisms themselves are universal, the combined one is also applicable to arbitrary systems, and it 
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can be deemed as a generalization. In Sec. E3, we first briefly review the two mechanisms respectively and then give 
the formulation of the competing mechanism. In Sec. Ill, we explain the differences between our mechanism and that 
adopted conventionally, taking ID Ising model as an example. In Sec. IV, as the chief task of this article, we apply 
it to kinetic Gaussian model and report the findings. In Sec. we summarize our study with some discussions. 



II. THE COMPETING MECHANISM 

First we briefly review the single-spin transition mechanism and spin-pair redistribution mechanism. 

A. Single-spin transition mechanism 

Glauber's single-spin flipping mechanism allows an Ising system to evolve with its spins flipping to their opposite. 
In single-spin transition mechanism [|ll|jl2| , a single spin Oi may change itself to any possible values, <7j, and the 
master equation is 



-P({a},£) = -53£[Wi("< - &t)P({<T},t) - Wi{&i -> ai)P({a jfH },&i,t)] , 



(1) 



The transition probability is in a normalized form determined by a heat Boltzmann factor, 



W i (p l — > <7j) = — exp 



<*j> 



. <3» = E exp 



-PHi Oi, E Oj 



(2) 



One can clearly see that this mechanism favors a lower Hamiltonian of the system. Based on the master equation 
Eqs.(0), one can prove that the time expectations of single-spin and correlation functions are 



dt 



(CT il (t)cr 4 2(t) • ■ -cr in (t)) 



-n{(Ti 1 a-i 2 ...<Ti n ) + ^2< E 



vi(^fe)=l 



P(Wht). 



(3) 



When n = 1, it is 



j t Qk (*) = -Qk (t) + E 
M 



^O-kWk (<Jk -> CTfe) 



P({a};t). 



(4) 



B. Spin-pair redistribution mechanism 

Kawasaki's spin-pair exchange mechanism allows an Ising system to evolve with its nearest neighbors exchanging 
their spin values. In spin-pair redistribution mechanism two neighboring spins, (Jjcri, may change to any possible 
values, <Jj(Ji, as long as their sum are conserved, and the master equation is 



d 



P({a},t) = E E l- W H ( a W -> tj*l)P({<r};-t) 

+Wji (a-j&i -> Ujoi) P {{0^,01^} ,Vj,&i;t)] 



(5) 



The redistribution probability is also in a normalized form determined by a heat Boltzmann factor, 



f 



Wji {ojOi -> bjdi) = —S aj+ai ^ j+a[ exp 



-fiHji (&j,ai, {<7 m } 



(6) 
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where the normalization factor Qji is 



Qjl = ^2 ^+<ri,»j+»i ex P —fftijl (^ji^h{ a m} 



Cj ,<Ti 



(Here it clearly favors a lower system Hamiltonian, but in the combined mechanism we shall change the sign before 
PUjl and make it turn to the opposite.) The time expectation of single-spin is 



dt 



w {a} 



,k-\-w 



P(W;t) 



(7) 



where d is the space dimensionality and means a summation taken over the nearest neighbors. 

C. The competing mechanism 

With the competing mechanisms, single-spin transition with probability p and spin-pair redistribution with proba- 
bility 1 — p, the master equation can be written as, 



dt 



P ({cr} , t) = pG me + (1 - p) K me , 



where the Glauber-type 



(8) 



(9) 



and the Kawasaki-type 



K me = E E {ajtn -> P (M ; t) 

+Wji -> <Jj(Ti) P ({a^j , cri^k} , a-j , <7/;i)] . 



(10) 



The Glauber-type mechanism is used to simulate the contact of the system with the external heat bath, and the 
transition probability is of the form given by Eqs.(j^). This mechanism favors lower energy of the system. The spin- 
pair redistribution mechanism is used to simulate the energy flux. The redistribution probability given above favors 
a lower system Hamiltonian, but what we need here is to the contrary. We can turn this tendency to its opposite if 
we change the sign before [3TLji in the redistribution probability. It has the following form, 



Wji {ojo-i -» 6- 3 ai) = —S aj+ai: a j+ai exp 



0Hji {a m } 



(11) 



This normalized form implies that the tendency is toward a higher system Hamiltonian. The word, competition, is in 
fact that between the two opposite directions either favored by one mechanism. 

As given above, we have already obtained the time expectation of single-spin with either mechanism respectively. 
One can prove that with competing mechanisms it will be, 



d 
dt 



(12) 



where the Glauber-type, 



Qk = -Qk (t) + Yl 



{a-} L <T fc 



^ CFkWk (dfe — > &k) 



(13) 
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and the Kawasaki-type, 



Qf = -2dq k (t) + J2J2 
w {a} 



,k-\-w 



P({a};t). 



(14) 



There are equations available for correlation functions with competing mechanisms, but in later studies we find this 
single-spin equation is enough to yield satisfying results. 



III. ON ISING MODEL 



As mentioned above, these years great efforts have been contributed to Ising model with competing dynamics. 
Because exact analytical treatment is too hard, most of the studies have been either approximation or Monte Carlo 
simulation. There are some differences between the conventional method and ours. In the following paragraph we 
present our considerations. 

The results one may expect directly depend on the expression of the transition (flipping, exchange, redistribution) 
probabilities. We think that there are two requirements: First, this probability should contain the Hamiltonian, 
and thus naturally favors either higher energy (Kawasaki-type) or lower energy (Glauber- type). Second, introducing 
temperature into it, we require that the transition be influenced by the heat noise. It is the first requirement that 
makes the two mechanisms compete and in all the studies it has been well adopted. However, due to the difficulties 
of actual practice, most of them used the simplified versions. In most of the studies, the temperature factor has not 
been introduced into the exchange probability, while in the flipping probability it has been only partly combined. 
Typically for ferromagnetic Ising model it has been set as: Glauber-type 

Wi = min{l,cxp(A Ei/KT)}, 

Kawasaki- type 



1, A Eij > 
0, A E i:i ^ 



In the new mechanism the transition and redistribution probabilities do not take the simplified versions. Besides 
some mathematical aspects such as normalization, their difference lies in the role of temperature. With the new 
mechanism there is greater influence of heat on the system. The fact, order for lower temperature and disorder for 
higher temperature, is not universally true in the phase diagrams obtained in earlier studies [Q-Q. However, we believe 
that this expectation will be unshakable if the system is governed by the new mechanism. In order to further study 
it, we applied our method to ID Ising model. It is well known that, due to the heat noise, there is only paramagnetic 
phase in II? Ising model. The analytical results we obtained confirms this conclusion, however one increases the 
energy flux, for all temperatures. (The details are in Appendix A, but we suggest you read it later for your easier 
understanding of our method.) ID Ising model governed by the conventional mechanism has been studied in Ref. 
H with approximation method and MC simulation. Further application on 2D Ising model is beyond the scope of 
this article. There is no better or worse, since the two mechanisms have different characteristics, but we think the 
comparison will be interesting and also feasible in practice. 



IV. ON THE KINETIC GAUSSIAN MODEL 



In this section we apply the new mechanism to the 3D kinetic Gaussian model and report our findings with the 
phase diagram. One and two-dimensional cases are quite similar. First we briefly review some basic properties of this 
model. 

Gaussian model, proposed by T. H. Berlin and M. Kac, at first in order to make Ising model more tractable, is an 
continuous-spin model. It has the same Hamiltonian as the Ising model (three dimensional), 

N 

— (3H = K ^ ^ (Tijk {(Ti+w.jk + Vij+w,k + &ij,k+w) , (15) 
i,j,k=l w 

where K — J/ksT and ^2 W means summation over nearest neighbors. Compared with the Ising model, it has two 
extensions: First, the spins a^k can take any real value between (— oo, +oo). Second, to prevent the spins from tending 
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to infinity, the probability of finding a given spin between <xyfc and Uij k + d&ijk is assumed to be the Gaussian-type 
distribution 

' b ( b 



f {(J ijk ) da l]k = y — exp ^- -a ijk I do- iifc , (16) 

where 6 is a distribution constant independent of temperature. Being an extension of Ising model, Gaussian model 
is quite different however. In the equilibrium case, on translational invariant lattices it is exactly solvable, and later 
as a starting point to study the unsolvable models it has also been investigated with mean field theory and the 
momentum-space renormalization-group method. 

As an example of the applications of single-spin transition and spin-pair redistribution mechanism, we have studied 
kinetic Gaussian model with both of them respectively. We summarize its dynamic properties as follows ]ll|-[T^| . The 
inherent dynamical competition of this model is that: the system tries to lower its Hamiltonian with the spins tending 
to infinity, while the Gaussian-type probability serves to restrict this tendency. Above the critical temperature, the 
prevailing heat noise permits only a disordered state, whereas below the critical point some kind of order will appear. 
Our study shows that, on translational-invariant lattices, the dynamic critical exponent z — 1/v — 2 is independent 
of space dimensionality and the governing dynamical mechanism. 

Now we turn to treat the 3D kinetic Gaussian model with the competing dynamics. ID and 2D systems can be 
treated in a similar way and they have qualitatively the same properties. In earlier studies we obtained the time 
expectation of single-spin. With the competing mechanism we can borrow these equations from Ref. |ll[] and [fl3|| . 

(1) With Glauber dynamics: 



d K 

jj.lijk W = Qijk = ~1ijk + ~^ + Qi+l,j,k + Qi,j-l,k + Qi,j+l,k + Qij.k-1 + 1i,j,k+l) ■ (17) 



(2) With Kawasaki dynamics: since we have changed the sign before system Hamiltonian in the redistribution 
probability, 

N 

+w,jk 0~ij-\-w,k 0~ij,k-{-w J j 

i,j,k=l w 



N 

=> (3H = (-K) 2^ a V k ( a i+w,jk + <Jij+ W: k + CTij^k+w) , 

i,j,k=l w 

all the expressions will remain the same if we switch K to —K. Here we will have to do the same, 
jj.Qijk (*) = Q%k = 2[b+ - Qijk) - (q%jk - qi-i,j,k)] 

+ [(<7i,j+l,fc — Qijk) — (Qijk — Qi,j-l,k)\ + lW,j,k+l — Qijk) — (qijk — QiJ,k-l)]} 

(-K) 

~*~2 \b + (-K)] ^ i ~ 1 '^' k ~ — <?i-l,j-l,fc) + (%Qi-iJ,k - q-ijk - ?i-2,i,fc) 

+ 2 (2qi + ij t k — qi+i,j + i,k — + (2Qi+l,j,k — qijk — Qi+2,j,k) 

+2 (2g , j j j_i ! )t — ffij-ijfc+i — 9i,j-i,fc-i) + (2<7ij-i,fe — qijk — Qi,j-2,k) 

+ 2 (2qij + i : k — Qi,j+l,k+l — qi,j + l,k-l) + (2qi,j+l,k — qijk — qi,j+2,k) 
+ 2 (2qi,j t k-l — <7i-l,j,fc-l — Qi+l,j,k-l) + (2qi,j,k-l ~~ qijk — qi,j,k-2) 

+2 (2q.ij t k+i — qi+i,j,k+i — q%-i,j,k+i) + (2qi,j,k+i — qtjk — qi,j,k+2)] ■ (18) 

This expression is in fact not as complex as it seems. Each term in a bracket is a second order derivative of q, and 
they will cancel each other if a summation is taken over all the spins. Thus X^jfc Qijk = 0- 
Then with the competing mechanisms, 

j t qijk(t)=pQ? jk + (l-p)Q? jk . (19) 
It is natural to first consider using the system Hamiltonian, 
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n 



J E 



(id) 



to characterize its behavior. If we find the Hamiltonian is decreasing, then the system is believed to be evolving to 
Ferromagnetic phase; if the Hamiltonian is stable around zero, the system must be in Paramagnetic phase; if the 
Hamiltonian is increasing, then the system is evolving towards Antiferromagnetic phase. Unfortunately, in Gaussian 
model although we can conveniently obtain the exact result for {cn<Jj) and the average value, Ylij ( a i a j) l-N 2 , we find 
it difficult to obtain an analytical result of the system Hamiltonian. In order to differentiate between these phases, 
instead we study the following two aspects. First, magnetization, 



M(t) 



(20) 



and second, 



ij k ij k 



\i+j+k 



Qijk (t) ■ 



(21) 



For pure ferromagnetic phase, if without any extra conditions, we expect nonzero magnetization, that is to say, 
M (t) ^ 0. At the same time, if the system is divided into two penetrating subsets, in one of which the spins have 
i + j + k being an odd number, and in the other one the spins have i + j + k being even, then one will find these 
two subsets are almost identical, and that leads to M' (t) = 0. So M (t) ^ and M' (t) = is the characteristic of 
ferromagnetic phase. For pure antiferromagnetic phase, we expect such a situation that the system can be divided into 
two penetrating opposing sublattices in the way mentioned above, one of positive spin and one of negative. Except 
the direction of the spins, these two sublattices are identical. This leads to M' (t) ^ but M (t) = 0. In disordered 
paramagnetic phase, we expect disorder of the whole system and both M (t) and M' (t) to be zero. Thus, if we get 
the evolution of M (t) and M' (t), we can decide in which phase the system is. 

With Eqs.(|Tg|), (|20| ) and ( fjl| ) we can write the time derivative of M (t) and M 1 (t) as, 



dM (t) 
dt 



1 \ - d . 

n E j t qk {t) 



(22) 



and 



dM' (t) 
dt 



W ? Jt q ' k {t) 



6K 
J, 



~P 1 + +6(1-P) 



QK — b 



b-K 



M'(t) 



(23) 



The details of Eqs.(|23j) are in Appendix B. Because we want K/b to be positive, the inequality 

(■> 6K \ or, .6K-b n 
- P 1 + - +6(1- p)t -^>0, 



becomes 



, -36 + Alp + V(1296 - 2808p+1561p 2 ) < — < 1, 
12p \ lb 



K 



(24) 



for < p < 1. 

Fig.l shows the phase diagram of 3D kinetic Gaussian model. It is divided into several regions by the following 
three curves, 



exp (—K/b) = e 



(25) 



exp (-K/b) = 1/e, 



(26) 
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and 



exp (—K/b) = exp 



^- (-36 + Hp + v/(1296- 2808p + 1561p 2 ) 



(27) 



Phase Diagrams for ID and 2D models are quite similar. 

We know that the critical point of 3D Gaussian model is K c = b/2d, where d is the space dimensionality. In the 
region above line (p5|), the temperature is higher than the critical value. Both M (t) and M' (t) are approaching 
zero exponentially, and this is identified as paramagnetic phase. The overwhelming heat noise permits no observable 
magnetization, and the Hamiltonian is static in equilibrium. When the temperature is below T c , the heat noise 
becomes secondary and some kind of order appears. 

For the region below line (^5|) and above curve d27j), as well as that below line (|26|), we have exponentially increasing 
M (t) while M' (t) is approaching zero, and this corresponds to ferromagnetic phase. In this phase, Glauber-type 
mechanism prevails and the energy is decreasing. There is observable homogeneous magnetization, the direction of 
which depends on that of the initial magnetization. 

In the region between curve (p6H and line ( p7f ) we have both exponentially increasing M (t) and M' (t) , and it can 
not be simply identified as ferromagnetic or antiferromagnetic. We give it a name, heterophase. In this region, the 
Kawasaki mechanism and the energy flux control the system, and the energy has a tendency to increase. The spin 
values are interesting. Our analytical results show that M (t) = M (0) e At and M' (t) — M'(0)e Bt . This means 
that, if the initial phase is ferromagnetic, M' (0) = 0, but M (0) ^ 0, later we will observe \M (t)\ increasing but 
M' (t) staying at zero. If initially the system is antiferromagnetic, M (0) = 0, but M' (0) ^ 0, later we will have 
\M' (t)\ increasing but M (t) staying at zero. If initially a disordered paramagnetic phase is given, M (0) = 0, and 
also M' (0) = 0, we will get both zero M (t) and M' (t). From this one may get confused, — what is the real picture? 
How is the energy sure to increase if the spin values are dependent on other conditions? Direct computer simulation 
reveals the key. For example, we apply the periodic boundary condition and initially set the system as, 

... + 1,+1,-1,+1,+1,-1,+1,+1,-1,+1, 

This leads to M (0) ^ 0, but M' (0) = 0. After 1 second, it becomes (we only give approximate values), 

. . . 0.7, 1.1, 0.2, 0.8, 0.6, 0.3, 0.6, 0.6, 0.2, 0.6, 
0.6, 0.2, 0.6, 0.6, 0.3, 0.6, 0.8, 0.2, 1.1, 0.7, . . . 



After 10 seconds, 



After 24 seconds, 



...9,13,5,17,3,17,3,15,6,11, 
11,6,15,3,17,3,17,5,13,9,... 



... - 249, 6941, -6735, 12089, -10041, 13228, -8902, 9922, -3755, 3435, 
3435, -3755, 9922, -8902, 13228, -10041, 12089, -6735, 6941, -249, . . . 

This shows the routine of self- organization the system chooses in this specific case, and it is a self-explanatory picture 
of how the energy manages to increase while M' (t) = 0. 

However, there is one exception. If at the beginning the system is set in a pure homogeneous ferromagnetic phase, 
we will find that the system remains in this homogeneous state. How can one explain this? Actually the system 
does have a tendency to increase its energy in this region, but it needs a hint to know exactly how. Here we cite 
the words of Tome and Oliveira M: It is "the result, in the case considered here, of a far from equilibrium process, 
namely, the continuous flux of energy into the system. Thus an instability of the usual (equilibrium) solutions leads 
the system toward states with spatial self-organized structure." The self-organization routine is closely related with 
a basic phenomenon of symmetry loss. A system under conditions that lead to paramagnetic phase is like a ball 
placed on the ground, for which all directions are identical. A system set in such a heterophase is like a ball placed 
on the top of another sphere. The upper ball has a tendency to fall down, but the direction depends on some kind of 
disturbance. If there is no disturbance at all, it will stay still. It will even have its energy increasing if the lower ball 
is being lifted up. Surely in reality we have never seen a ball stay still on the top of another, because disturbance can 
not be avoided. So discussion here is only to help qualitatively explain this strange phase and the mechanism hidden 
behind. 
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V. RESULTS AND SOME DISCUSSION 



In this article, we have given a systematic formulation of the new competing mechanism. The Glauber-type 
single-spin transition mechanism with probability p simulates the contact of the system with the heat bath and the 
Kawasaki- type spin-pair redistribution mechanism with probability 1—p simulates an external energy flux. These two 
mechanisms themselves are natural generalizations of Glauber's single-spin flipping mechanism and Kawasaki's spin- 
pair exchange mechanism. Thus, on the one hand, this mechanism is in principle applicable to arbitrary systems, while 
on the other hand, our formulation is able to contain a mechanism that just directly combines single-spin flipping and 
spin-pair exchange in their original forms (not simplified). Compared with the conventional one, the new mechanism 
does not assume the simplified versions. Their difference lies in the different role the system temperature plays. (As 
we have emphasized before, there is no better or worse, since they are different.) We applied the new mechanism 
to ID Ising model and used the analytical results to exemplify this difference. With the new mechanism, there is 
greater influence of temperature and the fact, order for lower temperature and disorder for higher temperature, will 
be universally true. 

In Sec. IV, we applied this mechanism to 3D kinetic Gaussian model. ID and 2D models can be treated following 
a similar way and have qualitatively the same properties. We used M (t) and M' (t) (their definition given in that 
section) to characterize the system. In the phase diagram of 31? model, we confirm the expectation, order for lower 
temperature and disorder for higher temperature. For temperature above T c , the system evolves to paramagnetic 
phase. For temperature lower than T c , we observe ferromagnetic phase and another heterophase (instead of the 
antiferromagnetic phase as guessed). This interesting heterophase is the result of the energy flux and self-organization. 
We have analyzed it in details in that section, and hope it will help make clear the self-organization process in phase 
transitions. With regard to Ising model, we also hope it will be a good reference when people are trying to solve the 
puzzling differences between the results yielded by approximation method and MC simulation. 
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APPENDIX A: CALCULATIONAL DETAILS ON THE KINETIC ISING MODEL 



We apply the new mechanism to ID Ising model. We derive the time expectation of single-spin using Eqs.(|l2]), 

d 



dt 



In this combined form, the Glauber- type term is as Eqs.(fL3l 



-Qk 



(*) + £ 



M 



d-kW k (<Tfe 



P({a};t) 



(Al) 



and the Kawasaki-type is as Eqs.(14) 



Q% = -2q k (t)+Y, 



P(W;t). 



(A2) 
(A3) 



First we calculate (Al) 



a k exp [Ka k (gfc-i + Ofc+i)] ~ Qfc exp [-Ka k (a k -i + Q-fc+i) 
exp [Kcr k (a k -i + cr fe+1 )] + exp [-Ka k (<r fe _i + a k +i)] 
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0, = -cr fe+ i 

tanh 2if, a^-i = &k+i = 1 
- tanh 2if, CTfc- 1 = crk+i = -1 

^ + CTfc + i)tanh2if. 



So 



^ffi^t (cTfe -» CTfc) 



P(W;i) 



(*) + ^ (fffe-i + <7fc+i)tanh2if. 



Second, we calculate flA2|), 



(A4) 



(A5) 



^ CTfcWfc,fc+l (cfcCTfc+i — » CTfcCTfc+l) 

g/c exp [-if (<Tfc-iCTfc + afc+igfc+2)] + gfc+i cxp [-if (ofc-iafc+i + CT fc CT fc+2 )] 

exp [-if ((T fc _l(7fe + CTfc + lCTfc +2 )] + CXp [-if (o-fc_:l.<7fc + l + o-fcO- fe+2 )] 
(<7fc + CTfc+l) /2, CTfc = CTfc+l 

tanh [-if (<7fc_i - cTfe +2 )] , cr fe = -<j k +i 

CTfc + CTfc+i (CTfc — CTfe + i) 



2 4 

(0fe_l - 1) (<7 fc+ 2 + 1) ((T fc _i + 1) (<7 k+2 - 1) 



-4 



tanh (-2if ) . 



CTfe+O-fc+l , 1 



(cfc_i — (Tk+2 ~ <Jk-\OkVk+\ + ^k^k+i^k+2) tanh (— 2if) 



And similarly, 



So 



^ &kWk,k-l (ofcOfc-l — > CTfcCTfc_l) 



' 7fc + 2 °' fc 1 + t (CTfe+i - o-fc-2 - o-fc+io-fccr/c-i + crfecrfe-icrfc-2) tanh (-2if ) . 



Qf = -2q k (t)+J2 



2_j °kWk,k+l (o-fcCTfc+l —> CTfcCTfc+l) 



+ ^kWk^k-l (CTfcCTfc_i — » CTfeCTfe_l) 



P(M;i) 



^ (<7fc-i + <?fc+i - 2g fe ) + ~ tanh (-2if ) (g fc+1 - + g fe _ a - <7 fc+2 ) 
+ - tanh (-2if) ((crfcCr fc+1 CT fc+2 ) - 2 ((J k +i^k^k-i) + ■ 



Thus, combining ( |A5| ) and ( |A6| ) one will get, 



(A6) 



~9k (*) + ^ + ftt+i) tanh2if 



1 



2 (ftfc-i + 3fe+i - 2%) + - tanh (-2if) (g fe +i - %_ 2 + g fe _i - qk+2) 



tanh(-2if) ((CTfcCT fe +io- fc+2 ) - 2 (CTfc+iCTfcCT fe _i) + {a k <J k -i(J k -2)) 



(A7) 
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In order to decide in which phase the system is, we have suggested two quantities in Sec. [V 



and 



We use their evolving tendency to characterize the system behavior. Obviously, (1) in a homogeneous ferromagnetic 
phase, we expect \M [t) \ /N — ► 1 and \M' (t) \ /N — > 0; (2) in an antiferromagnetic phase consisting of two penetrating 
and opposing sublattices, we will have \M' (t)\ /N — > 1 and |M (i)| /JV — > 0; (3) in disordered paramagnetic phase, 
both M' (i) /JV and M (i) /JV will approach zero. (A detailed analysis can be found in that section.) With Eqs.(A7) 
we can get, 

j t M (t) ee Ti qk W = ~P 0- - tanh 2tf) M (t) . 



Thus 



At the same time, 



M{t)=M (0) exp [~p (1 - tanh 2K) t] 



and 



(-) fffc (*) 



Vk W + 2(-9fe-i-^+i) tanh2 ^ 



+(l-p) 



<7k-l — 9fe+2 J 



+- tanh (-220 (- (<j' k <j' k+l( j' k+2 ) + 2 K^a^) - (p> k o> k _ 1 o> k _ 2 )) 



(A8) 



Thus 



M' (i) ee ^ — ^ (i) = [p (-1 - tanh 2K) + (1 - p) (-2 + tanh 22T )] M' (t) 

k 

= [-2 + p + (1 - 2p) tanh 22T] M ' (t) . 



M' (t) = M' (0) exp { [-2 + p + (1 - 2p) tanh 22T] i} . 



(A9) 



In ID Ising model, as given in Eqs.(AS) and (AS), M (t) and M' (t) are both approaching zero exponentially, and 
we can make the conclusion that, however one increases the energy flux, the system will stay in paramagnetic phase 
at arbitrary temperature. 



Kawasaki- type: 



APPENDIX B: CALCULATIONAL DETAILS OF EQS.(bd) 



dt q ' ijk ^ = 2{b-K) b ^ [(~ q id,k+i - Qijk) - (fltfk + 

+ [(-Qi+i,j,k ~ lijk) - (lijk + Vi-i,j,k)] + [(-1i,j+i,k - lijk) - (lijk + 1i,j-i,k)] } 

—K 

~^~2 (b — K ) ^ ( _ _ l'i-l,j+l,k ~ 1i-l,j-l,k) + \~^Qi-l,j,k ~ lijk ~ Qi-2,j,k) 
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+2 ( — 2q' i+1 - k — q' i+1 J + — q' i+1 j-i,k) + (~^i+l,j,k ~ l'ijk ~ 1i+2,j,k) 
+ 2 ( — 2g- J _ 1 fc — q'ij-lM+l — 9i,j-l,k-l) + ( — ^l'i,j-l,k — l'ijk — Qi,j-2,k) 
+2 (— 2q' i; j +l k — q'i^+i^+i — q'ij+i.k-l) + (~^Qi,j+l,k ~ Qijk ~ l'i,j+2,k) 
+2 (— 2q' i t j t k_i — q'i-ij,k-i ~ 9i+ij,fc+i) + ( — 2g- j^-i — 4ijk ~~ <7ij\fc-2) 
+2 (— 2g-^ fe+1 — <7i+i,j,jfe+i — Sj-i^jfc+i) + (—2<li,j,k+i ~ <2yft — lij.k+2)] 



Glauber- type: 



And then 



d K 

I i \ Iff f / / / f \ 

~^1ijk W = _ + ~g \ Qi-l,j,k ~ 1i+l,j,k ~ <lij + l,k ~ QiJ-l.k — 1iJ,k+l ~ 1i,j,k+l) 
dM' (t) 



dt 



k 

7< I 1 



6K 



(l-p)6 



6K -b 



b-k 



M' (t) 



(Bl) 



(B2) 
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Caption of figures 



Fig.l. The phase diagram of 3D Gaussian model with competing dynamics, Glauber- type with probability p, and 
Kawasaki- type with probability 1 — p. The system exhibits paramagnetic phase (Para), Ferromagnetic phase (Ferro), 
and heterophase (Hetero). 
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